Quantum dynamics of localized excitations in a symmetric trimer molecule 
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We study the time evolution of localized (local bond) excitations in a symmetric quantum trimer 
molecule. We relate the dynamical properties of localized excitations such as their spectral intensity 
and their temporal evolution (survival probability and tunneling of bosons) to their degree of overlap 
with quantum tunneling pair states. We report on the existence of degeneracy points in the trimer 
eigenvalue spectrum for specific values of parameters due to avoided crossings between tunneling 
pair states and additional states. The tunneling of localized excitations which overlap with these 
degenerate states is suppressed on all times. As a result local bond excitations may be strongly 
localized forever, similar to their classical counterparts. 

PACS numbers: 34.30.+h, 05.30.Jp, 03.75.Lm, 05.45.Mt 



I. INTRODUCTION 

The study of the classical and quantum dynamics of 
excitations in non-linear systems with few degrees of 
freedom has been used for decades to understand the 
processes of energy redistribution after an initial local 
bond excitation in polyatomic molecules Q, Q |^ ^ lEi IB| - 
Equally, interest in these issues evolved from the more 
mathematical perspective of nonlinear dynamics, local- 
ization of energy and solitons 0. This second path 
was boosted by the observation of discrete breathers 
(DB) - time-periodic and spatially localized excitations 
- in a hu ge v ariety of spatially discrete lattice systems 
H, 1^ Hol lllj. The flood of recent experimental obser- 
vations of DBs in various systems includes such differ- 
ent systems as bond excitations in molecules, lattice vi- 
brations and spin excitations in solids, electronic cur- 
rents in coupled Josephson junctions, light propagation 
in interacting optical waveguides, cantilever vibrations 
in micromechanical arrays, cold atom dynamics in Bose- 
Einstein condensates loaded on optical lattices, among 

others a^jj.allillilllllil[iamill2l. in 

a substantial part of these cases quantum dynamics of 
excitations is either unavoidable (molecules, solids) or 
reachable by corresponding parameter tuning (Josephson 
junctions, Bose-Einstein condensates). 

Progress in classical theory was achieved using a syn- 
ergy of analytical results and computational approaches. 
The computational aspect is vital because we deal gener- 
ically with non-integrable systems, which can not be 
completely solved analytically. Computational studies 
of classical systems of many interacting degrees of free- 
dom (say TV oscillators) are straightforward since we have 
to integrate 2N coupled first-order ordinary differential 
equations, so '--^ 10* is no obstacle to do even long- 
time simulations in order to study statistical properties. 
The quantum case is much less accessible by computa- 
tional studies. This is because in general each degree 
of freedom (e.g. an oscillator) is now embedded in an 
infinite-dimensional Hilbert space. Even after restrict- 
ing to only s states per oscillator, the dimension of the 



Hilbert space of the interacting system is now s^, making 
it nearly impossible to treat both large values of s and 
A^ - independent of whether we aim at integrating the 
time-dependent Schrocdinger equation or diagonalizing 
the corresponding Hamiltonian. However it is possible 
to treat small systems with A^ = 2, 3, which adds to the 
above mentioned studies of bond excitations in molecules, 
perspective cases of few coupled Josephson junctions, and 
Bose-Einstein condensates in optical traps with just a few 
wells [23, 01 • Remarkably, in the last case there is al- 
ready an experimental realization |25j |. 

Extensive studies of a dimer model N = 2 with addi- 
tional conservation of the number of excitations (bosons) 
have been accomplished jHESIIibg, 30, 31, 32, 33:|. 
The conservation of energy and boson number makes 
this system integrable. Due to the nonlinearity of the 
model the invariance under permutation of the two sites 
(bonds, spin flips etc) is not preventing from having clas- 
sical trajectories which are not invariant under permu- 
tations. These trajectories correspond to a majority of 
bosons (and thus energy) being concentrated on one of 
the sites. Quantum mechanics reinforces the symme- 
try of the eigenstates via dynamical tunneling in phase 
space (without obvious potential energy barriers being 
present) [3^ l35l | . The tunneling time is inversely propor- 
tional to the energy splitting of the corresponding tun- 
neling pairs of eigenstates. Notice that while most of the 
quantum computations concerned diagonalization of the 
Hamiltonian, a few results show consistency with numer- 
ical integration of the Schrocdinger equation 31., .3.7]. In 
|25l | the first experimental observation of non-linearity- 
induced localization of Bose-Einstein condensates in a 
double-well system was obtained, in agreement with re- 
sults discussed above. 

The extension of the dimer to a trimer A^ = 3 allows 
to study the fate of the tunneli ng p airs in the presence of 
nonintegrability [MIlllsllMISlii and in an effective 
presence of the fluctuation of the number of bosons (on 
the dimer). Diagonalization showed that tunneling pairs 
survive up to a critical strength of nonintegrability |37j , 
while the pair splittings showed characteristic resonances 
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due to interactions with other eigenstates |4ll |. 

Trimer models have been also extensively studied in or- 
der to describe spectral properties and then en ergy trans- 
fer in ABA molecules like water "4^, "i?, "S, "1^, "46] , which 
is connected with the appearance of quantum local modes 
(discrete breathers) . In these studies the presence of local 
modes was already identified as nearly degenerate eigen- 
states (tunneling pairs) in the eigenvalue spectra of the 
considered systems. 

In this work we study the time evolution of localized 
excitations in the trimer, and compare with the spectral 
properties of the system. We compute the eigenvalues 
and eigenstates of the quantum system and then the ex- 
pectation values of the number of bosons at every site 
on the trimer and the survival probability of different 
initial excitations as a function of time. We also com- 
pute the spectral intensity of the initial excitations to see 
how many eigenstates overlap are involved. That allows 
to draw conclusions about the correspondence between 
the time evolution of a localized initial quantum state 
(not an eigenstate) and the presence or absence of quan- 
tum breathers, i.e. dynamical tunneling eigenstates. We 
identify novel degeneracies in the trimer spectrum due 
to avoided crossings, and relate these events to unusual 
classical-like behaviour of quantum localized excitations. 

II. LOCAL BOND EXCITATIONS IN THE 
CLASSICAL CASE 

The classical trimer is described by the Hamiltonian 

M 

1 S 
H^Hd + -{Pi + Xi) + -{X1X3 + P1P3 

+X2X3 + P2P3), (1) 

Hd = ^iP? + Pi + Xl + Xl) + i[(p2 + xlf 

+ {Pi + Xlf] + j{XiX2 + P1P2), (2) 

where iJ^; is the dimer part. In all of this work we use 
dimensionless quantities. C is the coupling inside the 
dimer, and 6 is the coupling between site 3 and the dimer 
which also destroys the integrability of the system. Using 
the transformation ^f,; = {l/^/2){Xi + iPi) the Hamilto- 
nian becomes 

H = Hd + ***3 + + ^a^s + cc) (3) 
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FIG. 1: (Color online) Poincare map of the classical phase 
space flow of the trimer. The map condition is A13 = 0. The 
plotting plane is X — l^'ip and Y = |'I'2|^. The parameters 
are B = 40, C = 2, 5 = 1. 



is conserved, and hence the problem is effectively two- 
dimensional. Also the trimer (and the dimer) is invariant 
under permutation of sites 1 and 2. 

We are interested in the fate of localized excitations, 
where some energy is excited e.g. on site 1, and none 
on site 2 (inside the dimer). The third site may have 
some nonzero energy as well (like an environment). For 
different initial conditions 



*i(0) = yy + ^^, «'2(0) = 0, ^^{Q) = s^--u (5) 

we computed the time evolution of the quantities l^'ip — 
by numerically solving the equations of motion. In 
all computations we used B = 40, C = 2, and 5=1. We 
also generate a Poincare map (FigQ using the condition 
A13 = {^i = Aie"^'\ Ay = Lpi-ipj) and the plane X = 
I^E'ip, Y = |5'2p. We observe that for positive v the 
evolution is regular and not invariant under permutation, 
so most of the energy initially placed on site 1 stays there, 
with site 2 becoming only little excited. Negative values 
of v yield chaotic motion which is permutation invariant. 

This transition from localization to delocalization of 
energy is also nicely observed in the temporal evolution 
in FiglH Increasing v from negative to positive values 
the energy exchange between sites 1 and 2 of the dimer 
is stopped. 



+C7(*t*2 + cc), (4) 

and the equations of motion transform to i^',; = dH/d'i*. 
Note that the total norm B = ^J^fi -h 5'^*2 + *5*3 



III. LOCAL BOND EXCITATIONS IN THE 
QUANTUM TRIMER 

The quantum trimer is obtained after replacing the 
complex functions ^P* by the bosonic operators a and 
(rewriting = {l/2)(^'*^ + vE'vp*) previously to 
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FIG. 2: (Color online) Time evoluti on of {i = 1,2,3) 

for different initial states *i,3(0) = V20 ± u (vE'2(0) = 0): (a) 
ly = -6, (b) u = 0, (c) iy = 6. 

insure the invariance under exchange ^ <-!■ ^'*): 
3 

H = Hd + -ja\a3 + 5{a\a3 + a\az + c.c), (6) 

15 3 ' ' 1 

— + -{a[ai + ala2) + ^ [(alai)^ + (4a2)^] 

+C(ala2 + C.C.), (7) 

where we take h = 1. The boson number operator 
B = a\ai + a\a2 + <x\a^ commutes with the Hamilto- 
nian, so we may diagonaUze © in the basis of eigenfunc- 
tions of B, {|ni, n2, ns)}, where ni,n2,n^ respectively 
are the number of bosons at site 1, 2, and 3. There are 
(&+l)(6 + 2)/2 eigenstates in the subspace corresponding 
to a fixed value of the eigenvalue hoi B. Since the Hamil- 
tonian is invariant under permutation between sites 1 and 
2 we expanded the wave function in the basis of symmet- 
ric and antisymmetric eigenstates of B, {|ni, n2, n^)s,A\, 
where 

\ni,n2,n:i)s,A = ^(|ni, 712,713) ± |n2, 711,713)). (8) 



Then the initial state 1^*0) = |7io, 7710, /q) writes as 

l^'o) = -^(|no,7no,Zo)s' + |no,mo,/o)yi), 

= -L(|M/o)5 + |*o)a), (9) 

In this representation diagonalization of the Hamiltonian 
reduces to diagonalize two smaller matrices — symmetric 
and antisymmetric decompositions of H — whose eigen- 
values are Ejf''^\ with less computing cost than diago- 
nalization of the full Hamiltonian. All computations were 
done using this representation. 

We computed the time evolution of expectation val- 
ues of the number of bosons at every site on the trimer 
(7ii)(<) = (5't|7ii|^'t) and the survival probability Pt = 
|(5'o|^'t)P (see appendix for explicit expressions), start- 
ing with various boson number distributions among site 
1 and site 3 controlled by the number i/: j^I'o) = \b/2 + 

0, b/2 — v), with h — 40, C — 2, and 5 = 1. In computa- 
tions we dropped the two first terms of the Hamiltonian, 
which are diagonal and just shift the spectrum. 

Tunneling pairs and localization 

In Figl^l we show the time evolution of expectation 
values of the number of bosons in the trimer. When the 
initial excitation is mainly localized at the third site in 
the trimer there is a fast redistribution of bosons between 
the two sites in the non-linear dimer until the dimer sites 
are equally occupied ( FigUJ-a). As we place more bosons 
on the dimer (site 1) the tunneling time of the excitation 
increases rapidly until the time of computation becomes 
too short to observe slow tunneling. On these timescales 
we thus observe localization of bosons on one site in the 
dimer ( Fig|31-b and^-c), in analogy to the classical case. 
The reason for this behavior is the appearance of tun- 
neling pairs of symmetric and antisymmetric eigenstates 
with very close eigenenergies in comparison to the mean 
energy separation between eigenstates ( FigEJ. These 
pairs strongly overlap with the initial state, as observ- 
able from the spectral intensity = P in the inset of 
the FigH 

The results in Fig ^ show an enhancement of the sur- 
vival probability with increasing boson number at site 
1, which is consistent with the results discussed above. 
The dominant tunneling pairs in the spectral intensity 
(inset of FigQ give the main contribution to the time 
dependence of the survival probability [3^. 

Avoided crossings and degenerate eigenstates 

Energy levels exhibit avoided crossings when we vary 
the parameter 5 which regulates the strength of nonin- 
tegrability of the system [33, as shown in the Fig|Sl Of 
particular interest is the outcome of the interaction of a 
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FIG. 3: (Color online) Time evolution of expectation values 
of the number of bosons at each site of the trimer for different 
initial states |*o) = 120 + 1^, 0, 20 - i^): (a) p = -6, (b) u = 0, 
(c) z/ = 6. 



FIG. 4: (Color online) Survival probability of the initial state 
i*o) = |20 u, 0, 20-1/). (a) v = -6, (b) u = 0, (c) u = 
6. Inset: Spectral intensity of the initial state |^o). Filled 
circles — symmetric eigenstates, open circles — antisymmetric 
eigenstates. 



single eigenstate and a tunneling pair. The principal dif- 
ference between these states is that the member of a tun- 
neling pair has exponentially small weight in the dynam- 
ical barrier region, which is roughly defined by ni = n2 
in the ni — n2 plane. A single eigenstate will in gen- 
eral have much larger weight in this region. Since each 
eigenstate is either symmetric or antisymmetric, and a 
tunneling pair consists always of states with both sym- 
metries, the interaction with a third eigenstate will in 
general allow for an exact degeneracy of two states with 
different symmetries. While this is in principle possible 
for any two states of different symmetry, the exponen- 
tially small weight of the tunneling pair states in the 
dynamical barrier region makes a difference. Indeed, a 
linear combination of two states with large (not expo- 
nentially small) weight in a barrier region yields again, 
though an assymetric state, but one with large weight in 
the barrier region. Contrary, for the case of a tunneling 
pair and a single state, we may expect an asymmetric 
eigenstate which has much less weight in the barrier re- 
gion, leading to a much stronger localization of the state 




FIG. 5: (Color online) An enlargement of the trimer spectrum 
showing three particular pair-single state interactions (num- 
bered circles). Thin solid line — symmetric eigenstates; thick 
dashed lines — antisymmetric eigenstates. 
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FIG. 6: (Color online) An enlargement of the trimer spectrum 
around the avoided crossing 1 in figure |S] involving the anti- 
symmetric state A-268 (thick dashed line) and the symmetric 
states S-286 and S-285 (thin solid lines). 



similar to a classical one. 

We analyze three particular avoided crossings identi- 
fied by numbered circles in Fig|Sl by computing the en- 
ergy separation AE{5) between such a single state and 
a quantum breather tunneling pair. We identify three 
different situations. The first one shows that the energy 
levels intersect once in some degeneracy point (FigEJ. 
At some value of the parameter S the energy separation 
between one of the members of the tunneling pair and 
the single state vanishes. Tunneling is suppressed com- 
pletely, and then an asymmetric linear combination of 
the degenerate eigenstates will constitute a non-decaying 
localized state. This situation has been well described 
by perturbation theory |4l| , where effects of other eigen- 
states have been neglected. We computed the density 
p{ni,n2) = |(rii, ^2, nal^) p of the asymmetric eigenstate 

10) = Vl</'i'^^»/\/2, where |</.f ^) and are the 

degenerate eigenstates. The result is shown in the Fig|7| 
where we can see that there is only a partial localization 
of the excitation, since the wave function has visible con- 
tributions around the diagonal X = Y (ni =712). Note 
that in addition it also shows sizable contribution on the 
other side of the barrier (ni « 2, ^2 ~ 26 in FigEJ. In 
fact the expectation values of the number of bosons for 
this state are (m) = 14.99, (712) = 14.89, (713) = 10.12. 
Thus in terms of averages practically no localization oc- 
curs since (ni) « (772) despite the observable asymmetry 
in FigEI 

The other two cases appear as a consequence of the 
infiuence of other states in the spectrum. In one case 
the energy levels do not intersect at all ( FiglHJ, due to 
the presence of another avoided crossing located nearby. 
In the third case surprisingly we observe that the en- 
ergy levels intersect twice. The situation is shown in 
Figl^l and is sketched in the inset of the figure. Due to 
the interaction with other states of the system we ob- 
serve an intersection of the two states of the tunneling 
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FIG. 7: (Color online) Density of the asymmetric state j^) — 

(10285) + I^Ms))/"^ ^ ^ function of the number of bosons 
at sites 1 and 2 at the degeneracy point 5d- Here X is the 
number of bosons at site 1, and Y is the number of bosons at 
site 2. 
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FIG. 8: (Color online) An enlargement of the trimer spec- 
trum around the avoided crossing 2 in figure |H| involving the 
symmetric state S-287 (thin solid line) and the antisymmetric 
states A-269 and A-270 (thick dashed lines). 

pair at some distance from the actual avoided crossing 
with the third state. Consequently both states have ex- 
ponentially small weight in the barrier region, and we 
may expect a very strong localization. In Fig llOl we can 
see that the asymmetric quantum breather in the degen- 
eracy point S ~ 1.462 in Fig|^ (see arrow) is strongly 
localized and the tunneling is suppressed for all times. 
Note that in both cases two and three the order of the 
participating levels before and after the avoided crossing 
is not conserved, at variance to the first case we dis- 
cussed above and which was described also in reference 
[4l|. The abovementioned strong localization of this ex- 
act asymmetric eigenstate is reflected in the fact that 
the wave function has practically zero weight around 
the barrier region X = Y {ni = 712). Note that at 
variance with Fig[7| here the wave function has no siz- 
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FIG. 9: (Color online) An enlargement of the trimer spec- 
trum around the avoided crossing 3 in figure |S| involving the 
symmetric state S-284 (thin solid line) and the antisymmet- 
ric states A-267 and A-268 (thick dashed lines). The curves 
for antisymmetric eigenstates were generated from the data 
AE{S) by using Ea = Es-'yAE, with 7 = 1000, for better vi- 
sualization of intersection between eigenvalues. Inset: Sketch 
of the variation of eigenvalues participating in the avoided 
crossing. Thin solid line — symmetric eigenstate, thick dashed 
line — antisymmetric eigenstates. The arrows mark the ana- 
lyzed degeneracy point (see text). 
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FIG. 11: (Color online) Time evolution of expectation values 
of the number of bosons at every site on the trimer at the 
degeneracy points for the cases shown in (a) Fig|71 (b) FiglTUl 
The initial state is |*o) = |26, 2, 12). 
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FIG. 10: (Color online) Density of the asymmetric state 10) = 

(I "^284) + \ 4'268)) / *s ^ function of the number of bosons at 
sites 1 and 2 at the degeneracy point Sd- Here X is the number 
of bosons at site 1, and Y is the number of bosons at site 2. 



able contribution on the other side of the barrier as 
well. This state is thus very close to its classical dis- 
crete breather counterpart (Fig^. Indeed, for this state 
(ni) = 25.62, (na) = 2.38, (^3) = 12.00. Consequently 
we find a very strong localization for the expectation val- 
ues, in addition to the observed asymmetry in Fig llOl 

It is interesting to test whether initial states with 
some distribution of bosons at every site of the trimer 
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FIG. 12: (Color online) Survival probability of the initial state 
1^0) = |26, 2, 12) at the degeneracy points for the cases shown 
in (a) FigUl (b) Fig llOl Inset: Spectral intensity of the ini- 
tial state 1^*0). Filled circles — symmetric eigenstates, open 
circles — antisymmetric eigenstates. 
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(l^'o) = |no, Too, /q)) can significantly overlap with the 
above described asymmetric eigenstates. This distribu- 
tion is given by the maxima in the density for every asym- 
metric eigenstate (around ni — 26, and 712 ~ 2). In 
figures ^2 ^^nd ^] we show the time evolution of the ex- 
pectation value of the number of bosons at every site on 
the trimer and the survival probability of such an initial 
excitation. A detailed analysis of the spectral intensity 
of the initial state |*o) = |26, 2, 12) (inset in the FigUH) 
shows that this initial excitation overlaps strongly with 
the degenerate eigenstates corresponding to the strong lo- 
calization shown in Fig llOl It implies that this excitation 
( Fig llU b) will never distribute its quanta evenly over 
both sites of the dimer. For the case shown in FigllUa 
the initial excitation has a smaller overlap with the de- 
generate eigenstates which gives the partial localization 
shown in Fig|7| Since the overlap is not zero the excita- 
tion will also stay localized in the sense that the crossing 
of curves corresponding to (rii) and (712) as in Fig|31-a and 
will never occur. Note that despite the difference be- 
tween the analyzed cases one and three, the evolution of 
the expectation values and the survival probabilities do 
not differ drastically. It needs more sensitive details in 
the preparation of an initial state to observe a practically 
total localization of bosons on one of the dimer sites for 
case three as compared to case one. 



IV. CONCLUSIONS 

In this work we observed how spectral properties of 
the Hamiltonian are reflected in the time evolution of 
different localized excitations in a trimer molecule model 
by monitoring the spectrum, the time evolution of ex- 
pectation values of the number of bosons at every site on 
the trimer and survival probabilities of different localized 
excitations. 

The tunneling pair splitting determines the lifetime of 
localized excitations. The survival probability and the 
time evolution of the expectation values of the number 
of bosons are clear indicators for a localized excitation 
being close or far from a quantum breather tunneling 
pair, while the spectral intensity of localized excitations 
is typically broad and does not show the peculiarities of 
the tunneling dynamics. Probing the time evolution of 
initially localized excited states thus allows to conclude 
about the presence or absence of tunneling pair eigen- 
states. 

We report on the existence of degenerate levels in the 
spectrum due to the presence of both avoided crossings 
and tunneling pairs. In these degenerate points tunnel- 
ing is suppressed for all times. While in general the 
asymmetric exact eigenstates will have quite large weight 
in the dynamical barrier region (contributed by the sin- 
gle level), we observe specific parameter cases where the 
weight is very small and the corresponding asymmetric 
eigenstate very strongly localized. Full or partial localiza- 
tion of bosons appears for all time scales for some spe- 



cific states and some specific values of the parameters. 
This effect could be studied in experimental situations 
of Bose-Einstein condensates in few traps which weakly 
interact, as well as in systems of few coupled Josephson 
junctions which operate in the quantum regime. Tun- 
ing experimental control parameters will allow to lock 
localized excitations for specific values and both prevent 
the excitation from tunneling, as well as allowing for a 
fine tuning of the tunneling frequency from a small value 
down to zero in the vicinity of these specific control pa- 
rameter values. 



Acknowledgments 



We thank G. Kalosakas, S. Keshavamurthy, M. 
Johansson, and L. Schulman for useful discussions. 



APPENDIX: EXPECTATION VALUES AND 
SURVIVAL PROBABILITY 



Expanding the wave function in the basis of symmetric 
and antisymmetric eigenstates of the Hamiltonian 



^^.(^)^..^<^.,|^.(^)^^ (A.l) 



where ^^^^'^^ = {c^^^^'^^^^o) s,A and rnf,n^U = 
{4>^\ni,n2,n^) S.A, the expectation value of the number 
of bosons at site i writes as 



{nm = {n\'^){t) + (n^OW + K^>W, (A.2) 
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where 



a{S,A)jlJ-'{S,A)i{E^ 



xF 





{S,A) 



(S,A)_^(S,A)^j 



(A.3) 

(A.4) 



{ni}s,A 



xG 



iS,A) 



{"i}s,^ 



^ Vni,n2,ra3 



(A.9) 



(A.10) 
(A.11) 



(nr)(i) 



(A.6) 
(A.7) 



— 712) n2,ri3' (-^•^) where bars mean complex conjugation. 
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